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Abstract 

A major goal of the BepiColombo mission to Mercury is the determination 
of the structure and state of Mercury's interior. Here the BepiColombo rotation 
experiment has been simulated in order to assess the ability to attain the mission 
goals and to help lay out a series of constraints on the experiment's possible 
progress. In the rotation experiment pairs of images of identical surface regions 
taken at different epochs are used to retrieve information on Mercury's rotation 
and orientation. The idea is that from observations of the same patch of Mercury's 
surface at two different solar longitudes of Mercury the orientation of Mercury 
can be determined, and therefore also the obliquity and rotation variations with 
respect to the uniform rotation. 

The estimation of the libration amplitude and obliquity through pattern match- 
ing of observed surface landmarks is challenging. The main problem arises from 
the difficulty to observe the same landmark on the planetary surface repeatedly 
over the MPO mission lifetime, due to the combination of Mercury's 3:2 spin-orbit 
resonance, the absence of a drift of the MPO polar orbital plane and the need to 
combine data from different instruments with their own measurement restrictions. 

By assuming that Mercury occupies a Cassini state and that the spacecraft 
operates nominally we show that under worst case assumptions the annual lib- 
ration amplitude and obliquity can be measured with a precision of respectively 
1.4 arcseconds (as) and 1.0 as over the nominal BepiColombo MPO lifetime with 
about 25 landmarks for rather stringent illumination restrictions. The outcome 
of the experiment cannot be easily improved by simply relaxing the observational 
constraints, or increasing the data volume. 



1 Introduction 

Mercury is the innermost and smallest planet in the Solar System, orbiting the Sun 
once every 87.969 days. Comparatively little is known about Mercury. Ground-based 
telescopes reveal only an illuminated crescent with limited detail and are restricted to 
the moments of maximal solar elongation of the planet as a result of its proximity to 
the Sun. Mercury is also a difficult target for in-situ observations due to the hostile 
radiation environment and, because of its deep position in the gravitational well of the 
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Sun, the difficulty of putting a spacecraft into orbit around it. The only in-situ meas- 
urements available up to the present day come from only two spaceprobes, Mariner 
10 and MESSENGER. Mariner 10 mapped about 45% of the planet's surface in three 
fly-by 's in 1974 a r id 197 5, and discovered a strong and unexpectedly dipolar magnetic 
field (|Ness et al.1 . Il975l ). NASA's MESSENGER spacecraft already performed three 



flyby's of Mercury (2008-2009) bringing the surface mapped to an astonishing 90%. 

Mercury is currently the subject of renewed attention through the NASA MES- 
SENGER probe but also through the upcoming joint ESA/JAXA BepiColombo mis- 
sion. The latter consists of two probes (MPO and MMO) that will orbit Mercury and 
are scheduled to be la unched in 2014 and inserted into different orbits in July 2020 
(|Benkhoff et al.l . l2010l ^. The Mercury Planetary Orbiter (MPO) is situated in a low 
altitude orbit optimised for planetary survey is mainly dedicated to study the surface 
and the internal structure of Mercury. BepiColombo's scientific objectives focus on a 
global characterization of Mercury through the investigation of its interior structure 
and composition; surface; exosphere composition and dynamics; structure, dynam- 
i cs and origin o f itsmagnetized envelope; and testing of Einstein's general relativity 
( BenkhofF et al. . 2O10l ). One spacecraft, the Mercury Planetary Orbiter (MPO), is led 



by ESA and its pay load comprises eleven experiments and instrument suites. The 
MPO will focus on a global characterization of Mercury through the investigation of 
its interior, surface, exosphere and magnetosphere. In addition, it will test Einstein's 
theory of general relativity. The second spacecraft, the Mercury Magnetosphere Or- 
biter (MMO), is led by JAXA and will carry five experiments or instrument suites to 
study the environment around the planet including the planet's exosphere and mag- 
netosphere, and their interaction processes with the solar wind and the planet itself. 
(|Benkhoff et al.l . hoi& i 



Here we are mostly concerned with Mercury's internal structure, which is the 
most puzzling among the terrestrial planets and for which the Bep iColombo space 



mission will play an important role in advancing our understanding (jBenkhoff et al. 



2010l ^. Both missions are under way to play th eir distinct ro l e in th e roadmap to our 



understanding of Mercury as explained by e.g. iBaloeh et al.l (j2007l ) 



Mercury has a slow rotation period of 58,68 days and revolves around the Sun in 
87.96 days, corresponding to a 3:2 spin-orbit resonance. It therefore completes three 
rotations about its rotation axis for every two orbits. This spin-orbit resonance is 
unique in the solar system as all other occurrences of spin-orbit resonance in the solar 
system involve a planetary satellite and most happen to be in 1:1 resonance around 
their central body. The orbit of Mercury has the highest eccentricity of all the Solar 
System planets, and it has the smallest obliquity (axial tilt with respect to the orbital 
plane) . 

Since Mercury has an equatorial asymmetry, the Sun exerts a gravitational torque 
on it. The torque causes a periodic variation of the planet's rotation with a period of 
88 days, which is called the annual forced libration. The amplitude of this libration is 
linked to the state of Mercury's core and so an elegant way to obtain information on 
the deep interior not directly accessible to us is through its effect on a comparatively 
easily observable quantity. In the presence of a liquid core, the amplitude of the annual 
forced libration is about two times larger or more than in the case of a completely solid 
planetary interior because the liquid core does (almost) not participate in the libration 
and the moment of iner t ia of t he mantle is only half or less of the total moment of in- 
ertia of Mercury (jPeald . Il974l : iRivoldini et~all . l2009l l. The magnetic field observed by 



Mariner 10, if caused by dynamo action, is evidence for a liquid outer co re, which very 



likely contains a solid inner core according to thermal evolution models (jBreuer et al. 
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20071') ■ The r e cent d etection of an annual libration amplitude of 35.8 it 2.0arcsec by 
Margot et all (|2007l ') through Earth-based radar measurements, seems to confirm the 
existence of a liquid core. Besides the state of the core, libration observations can 
also be used to constrain the core size and composition. The large core of Mercury, 
as indicated by Mercury's high mean density, is thought to consist of iron and lim- 
ited amoui its of lighter elements, of which sulphur is thought to be the most likely 
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Van Hoolst and Jacobs. 2003; Schubert et al. 



candidatej Harder an d Schubertl. 

19881 : ISpohn et alT EoOl: Rivol dini et al.l . 120091 '). If Mercury has formed close to its 



curr ent position, i.e. close to the Sun, its c ore sulphur concentration is probably very 
low (jLewisl . Il988l : ISohl and Schubert! l2007l ^ . In the event that Mercury has formed in 
the same zones as Venus, Earth and Mars, its light element concentration could be a 
lot higher and close to that of the other three terrestrial planets the Earth has a core 
sulphur concentration of about lOwt.% , and for Mars a light element concentration 
of 14 wt.% is often considered ( Van Hoolst . 20071 ) (wt.% = weight percentage). 

The light element concentration of the core is related to the size range of the core 
and thus to the principal mantle moment of inertia. For larger cores and smaller mantle 
moment of inertia, the libration amplitude increases since the libration is essentially a 
movement of the mantle. The angular libration amplitude at the equator is about 21 
arcsec in the case of a completely solid core, and between 36 and 78 a rcsec for models 
with a liquid core with between 0.01 wt% and 8.6 wt% of light element (jRivoldini et al. 
20091 ^. 



Also the obliquity c an be used to constrain the interior of Mercury, as can be easily 
seen from Peale's idea ( Peale . 19761 ) to combine the measurements of the amplitude 
of the 88-day annual libration amplitude 7, the obliquity rj and the second degree 
gravitational field coefficients C20 and C22 into the equation: 
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where and i?^ are the mass and radius of Mercury, A < B < C are the planet's 
principal moment of inertia. C is aligned along the axis of planetary rotation (z- 
axis) while the x and y-axis correspond to the A and B principal moment of inertia, 
respectively in the equatorial plane. Cm is the polar principal moment of inertia of 
the solid silicate part of the planet (mantle -|- crust) alone and is aligned with respect 
to the z-axis. The first factor in the left hand side of eq.([T]) follows from the amplitude 
of the annual libration which is proportional to {B — A) as a result of the Solar torque 
and inversely proportional to Cm- The second factor is four times the gravitational 
coefficient C22: 
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and the third factor follow s from the Cassini state 1 equilibrium obliquity satisfying 
the equation (Peale, 198ll ): 
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where ijc is the obliquity of Cassini state 1, / is the inclination of Mercury's orbital 
plane with respect to the Laplace plane and /xlp is the precession rate of the orbital 
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plane about the Laplace plane normal (see also Yseboodt and Mareot . 20061 ) . As the 
orbital characteristics of Mercury are accurately known eq.([T]) shows that the obser- 
vation of the obliquity, libration, and degree-two gravitational coefficients is sufficient 
to determine Cm/C, the ratio of the mantle to the polar moment of inertia. This 
reasoning is only valid if Mercury is in or c lose to Cassini state 1. The spin axis ori- 
entation determined bv lMareot et alJ (120071) lie s in th e plane defined by the orbit pole 
and the Laplace pole of 
Recent observations by 



Yseboodt and MargotI ( 20061 ) as required for a Cassini state. 



Margot et al.l (j2007l ) confirm the validity of this hypothesis 



with a measured obliquity val ue of about 2.11 arcmin, which lies precisely within the 
theoretically predicted range (Peak, 19881 : Peale et al. . 2002). 

The main goal of the BepiColombo rotation experiment is to obtain insight into 
Mercury's deep interior from observations of variations of Mercury's uniform rotation 
and precise measurements of its orientation in space. In particular this experiment 
uses observations of the 88-day forced li bration, the obliquity, a.nd th e degree-two 
coefficients of the gravity field of Mercury (Peak , 19761 : Peale et al. . 2002). The gravity 
coefficients will be determined fror n the radio scie nce experiment MORE (Mercury 
Orbiter Radio-science Experiment) ( less et al. . 20091 ) while the obliquity and librations 
will be determined from the comparison of observed target surface positions at different 
times. 

Using the on-board star tracker to precisely determine the MPO's orientation and 
using the position determined through the radio-science experiment MORE, the loca- 
tion in inertial space of each image taken from the planetary surface by the SIMBIO- 
SYS HRIC (Spectrometer and Imagers for MPO BepiColombo Integrated Observ atory 



System High Re solution Imaging Channel. dThe SIMBIO-SYS Team et al.l . boid 'l) 



be determined (jless et al.l . l2009l : iMilani et al.l . l200ll ) , and changes in landmark posi 



can 



tions can be measured. 

The realistic estimation of the attainable accuracy of both the obliquity and lib- 
ration amplitude determination in the fran ie of the Bep i Colorn bo n iission is the main 
focus of this work. Previous simulations by Milani et al. ( 2001 ) and Jehn et al. ( 20041 ) 
used the best knowledge of the BepiColombo mission available at that time. Today 
many mission parameters potentially influencing the rotation experiment have changed 
and new constraints arisen. Additionally we also take into account the main gravita- 
tional perturbations the MPO will be subjected to. 

Such an analysis requires the use of an accurate model of the rotation of Mer- 
cury, which takes into account other longitudinal librations besides the main 88-day 
libration due to planetary perturbations on Mercury's orbit. As the long-period lon- 
gitudinal libratio n amplitudes depend on the same inte rnal structure quantities as the 
annual libration (jPufev et al.1 . 12OO8I : IPeale et al.l . I2OO9I ') we study the possibility of an 
independent determination of each amplitude in order to better constrain Mercury's 
interior structure. 

The paper is organised in the following way: In the following section the different 
models used to perform this study will be described. First the Mercury rotation model 
used (sec. 12. ip . then the MPO orbit model (sec. 12. 2p . the observable and the method 
used to simulate the rotation experiment (sec. 12. 3p and the considered error models 
(sec. 12. 4p . We next discuss the possibility to perform repeated surface observations 
from the MPO in sec. [3l In section H] we discuss the achievable precision levels on 
the annual and planetary librations as well as the obliquity. Finally some concluding 
remarks and perspectives are laid out in section [5l 
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2 Method 



In the frame of this work, we create a simulated dataset consisting of the central 
position of potential surface images as a functi on of time as seen b y the BepiColombo 
MPO. We consider an MPO orbit as defined bv lCarcia et al.l ( 200?! ) (see table[2]) taking 



into account the perturbative effects of the low-level gravitational field coefficients C20, 
C22 and C30 (see section [2T2|) . By projecting the spacecraft orbit on the surface of 
the planet, the distribution of potential images on the surface is then determined 
(see section [23]) . Each potential image is subjected to a positioning error due to the 
uncertainties on the satellite position and orientation. The image-pairing process is 
represented by an additional error depending on the photographic resolution of the 
least well-resolved image and strict limits on the allowed surface solar illumination 
angle, (see section [2^ The simulated dataset is then used to estimate the accuracy of 
the reconstruction of the orientation and rotational motion of the planet, as a function 
of the quantity of measurements made, the number of different targets considered and 
their locations on the surface of the planet in sections 3 and 4. 

2.1 Mercury Rotation Model used 



Peale et al.l ([200^) have shown that the differential equation governing Mercury's lib- 



ration 7 can be expressed as 



(PM 



-aY,G2Qq[e) sin [27 + (1 - g)M] 



(4) 



where M is the planetary mean anomaly and the G20q{e) are eccentricity functions of 
the type Gipq{e) as defined bv iKaula (Il966l . (chap. 3)) and a — - 
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The largest libration has a period equal to the period of the orbital motion of the 
planet around the Sun (87.969 days = 1 Mercury year) and can be obtained from 
eq-dH) by keeping only those terms proportional to cosM or sinM. This condition 
imposes that g = or g = 2 in eq.dH) yielding: 



7 = a (G200 - ^202) • (5) 

In general Mercury's forced libration can be expressed as a series of sine terms repres- 
enting the harmonics of the annual frequency. In order to describe Mercury's libration 
accurately enough for our simulation purposes, we need to take into account all har- 
monics with an amplitude larger than the expected precision of the BepiColombo 
libration experiment. It suffices to take the semi-annual term into account: 



7(t) = 7sinM(t) +744sin2M(t) , (6) 

where 744 is the semi-annual amplitude. Incidentally the semi-annual libration amp- 
litude does not add an unknown to our system, as the amplitudes of the higher har- 
monics of the annual libration are related through a fixed coefficient to the annual 
libration amplitude. These coefficients, which we call K for the first harmonic, only 
depend on the eccentricity of Mercury's orbit through the eccentricity functions. In 
the case of the semi-annual (44-day) libration we have: 

_ 744 G20-I — 6*203 ^y-j 

7 4(^200 - G202) ' 
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where the eccentricity function Gi^p^q{e) only depends on the eccentricity of the Mer- 
cury's orbit. For Mercury's current eccentricity value of e = 0.2056 this leads to a 
value oi K = —0.105455. The semi-annual amplitude of Mercury's libration is of the 
order of 10% of the annual amplitude, and we can write: 



'y(t) = -f{smM{t) + Ksm2M{t)) 



(8) 



Through the same reasoning we obtain that all amplitudes of shorter period harmonics 
of the annual forced libration will be smaller than 1 arcsecond. As we will show in 
section m the accuracy needed on the libration amplitude to fulfil the minimum mission 
requirements corresponds to ~ 3 arcseconds. Therefore there is no need to model the 
forced libration to a higher precision in the frame of this work. 

Long period terms (from which M is absent) can be obtained by averaging eq.(j3]) 
with respect to M over a whole revolution. We only have terms from which M is 
absent if q = 1. Since torques are small the angle 7 is considered very small and 7 
can be considered constant over one orbital period. Therefore, time-dependent forcing 
terms proportional to sin 27 can be neglected in the right hand side of eq.Q, reducing 
the equation to the harmonic oscillator equation. The resulting solution corresponds 
to the period of Mercury's long-term libration (or free libration, npree): 



_i 27r 

Hfree = (2aG'2oi) 2 — 
n 



(9) 



( Jehn et al.l . 2004). For a value of a of 3.06 x 10' 



obtained by Margot et al. ( 20071 )) and the expression of G201 = — ^e'^ this long- 



(based on the libration estimate 
7^ 123 .=3 
16 



term period corresponds to ~ 12 Earth-years. Due to different dissipation processes 
inside the planet this librati on should ha ve damped out and no observable amplitude 

Peald (|2005l l has shown that the damping time is < 



should have remained, since 



10^ 



years for all types of known excitations. This reasoning is valid except if a hitherto 
unknown process excites this libration. 

Planetary induced lihrations are caused through the same interaction mechanism 
as forced librations. The periodic torque arising from the interaction of a certain Solar 
System body with Mercury will cause a perturbation of Mercury's orbital motion, al- 
tering its position relative to the Sun and thus the solar torque acting on its equatorial 
bulge. This in turn will modulate the planet's libration motion and can be modelled 
as an additional libration of the same period as the perturbing forcing. The most 
important libration is caused by Jupiter and has a period of 11.86 years, coinciding 
with the time it takes for Jupiter to orbit the Sun once. The main libration caused by 
Venus has a period of 5.66 years correspondi ng to a 2:5 resonance between Venus' and 



Me rcury's solar longit udes (orbital periods) ( Peale et al. . 2007 : Dufev et al. . 20081 ) 



Peale et al.l (|2009l l have shown that for the accepted value of a = 3.06 x 10 * 
Mercury's libration will be dominated by two periods, the 88-day forced libration 
period and a long period libration forced at Jupiter's orbital period. The excitation 
mechanism at the basis for the large Jupiter induced libration amplitude is to be found 
in the proximity with the planet's free libration period. The independent measurement 



B-A 



because 



of its amplitude will pose an additional and independent constraint on -tt 
its amplitude is highly sensitive to the free libration period. 

To perform the Peale experiment Mercury needs to occupy Cassini state 1, as 
eq.(|3l) is based on this hypothesis. The recent Earth-based radar observations by 



Margot et al.l (|2007l ) have yielded a fairly accurate measurement of Mercury's obliquity, 
and indicate that this is the case. 
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If we define 70 as the initial angle 7 of the prime meridian at time to = J2000.0, 
and make use of the fact that the uniform rotation rate of the planet is equal to |M, 
the complete model of Mercury's rotation can be expressed as: 



e (t) = 7 + Mt + 7(t) + 7^ (t) + 7$ (t) , (10) 

where 77,. (t) and 75 (t) represent the main Jupiter and Venus induced libration re- 
spectively. as defined above is expressed in Mercury's non-rotating equatorial frame 
(Hermean mean equator of J2000.0). The detailed form of the planetary librations can 
be expressed as follows: 



79(0 = 7$ sin(2Ao - SA^ + (^5) 



(11) 



where Aj and ipi represent the planetary longitude and libration phase of planet i and 
Ao is Mercury's longitude. 



2.2 MPO orbit model 

In order to precisely model the trajectory of the BepiColombo MPO spacecraft around 
the planet we have to take into account several deviations from the perfect motion. 
Since the planet cannot be considered as a spherically symmetric body, for example 
because of the equatorial bulge, various orbital perturbations arise. Hence, we include 
effects on the spacecraft orbit of the low degree gravitational field harmonics C20, C22, 
C30. As the value of C30 has not yet been detern iined for Mercury w e use the following 



2003). We have verified 



estimate: C30 = O.IC20, as is also done by ESA (jGarcia et al. 
that using a value of opposite sign for C30 does not impact the results of this study in 
a significant way. 

The value used for the equato rial fiatteniiig of t he gravitational potential is the 
MESSENGER fivbv estimate from lSmith et al.1 (f ioinl 'l . Coo = 0.81^0.08 x 10"^ The 
valu e of the polar flat t ening of the gravitational potential, C20 = 6.0ib2.0 10~^ is taken 
from lAnderson et al.l (|l987l ) and based on the Mariner 10 mission flyby's. The value of 
C20 recovered from the MESSENGER flybys, is considerably lower in magnitude than 
the value obtained from Mariner 10 tracking. However, it is not as well constrained 
as the equatorial ellipticity because the fl y by trajectories wer e nearly in the planet's 
equatorial plane. Additionally the value by Smith et al. ( 2O10l ) is inconsistent with the 
Cassini state assumption. In order for Cgn to be fully consistent with the Cassini state 
equation when using the Smith et al. ( 20ld ) value for C22, its value should be taken 



Body 


Symbol 


Period 


Forcing argument 


Normalized Amplitudes 










Peale et al. (2009) 


Dufev et al. (2009) 


Sun 





43.98466 d 


2Ao 


0.10223 


0.11150 


Sun 





87.96935 d 


Ao 


1.00000 


1.00000 


Venus 


? 


5.66316 y 


2Ao — 5Xv 


0.09697 


0.10691 


Jupiter 




11.86295 y 


Aj 


1.09339 


1.2193 



Table 1: A summary of the main libration amplitudes and associated periods, normal- 
ized relative to the annual libration amplitude. No other forced libration is modelled 
in this study 
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inside the following interval —5.6 x 10~^and — 4.0 x 10~^, w h ich co rresponds to the 
lower part of the confidence interval given bv lAnderson et al.l (jl987l ). We choose the 
value for C20 to be —4.9 x 10~^ and for C22 to be 0.81 x 10~^, in our simulations. 



We calculate the motion of the spacecraft numerically integrating the Lagrang e 
equations which govern the variation of the osculating orbital elements ( Kaula . 19661 ). 
This process assumes a trajectory influenced by Mercury's low level gravitational field 
only (020,^22,(^30) and does not take into account for potential orbital manoeuvres or 
spacecraft wheel offloading. We have tested that 6*20,^22 pairs at the boundaries of 
the expected range do not affect the results. 



2.3 Observable 

In order to perform the simulation of this experiment we need to express the position of 
any given image of the planetary surface in an inertial frame. This will then enable to 
determine the rotation parameters of the planet from comparing the image positions 
using landmarks in them. In this study this approach is simplified by comparing 
surface point positions directly. 

A point on the planetary surface can be expressed in a Mercury body-fixed frame 
in spherical coordinates (r, <j), A), where r is the planetary radius, (j) is the latitude and 
A is the longitude of the surface point. The spherical coordinates are defined with 
respect to the rotation axis z, the x-axis pointing in the direction of the planetary 
prime meridian at J2000 and the y-axis in the equatorial plane perpendicular to the x 
and z-axis. We here assume Mercury's shape to be spherical, to simplify our surface 
position calculations. This simplification amounts to the hypothesis that the shape of 
the planet is well known and does not represent an error source for the experiment. 
We express the position of a fixed point on the planetary surface in an inertial frame 
by performing a series of rotations. First a rotation from the rotating body-fixed frame 
to a Mercury non-rotating frame is performed. We then perform a transition from an 
equator-based reference frame to an orbital-plane-based reference frame through the 
application of 3 rotations with respect to the axes z-x-z, in order to take the planetary 
obliquity r/ into account. The first rotation about the 2-axis aligns the equatorial 
frame's x-axis with Mercury's vernal point at J2000 . The second rotation about the 
X-axis, then aligns the equatorial plane with the orbital plane, with an angle equal to 
the obliquity. The third rotation is an inverse of the first in order to align the x-axis 
with Mercury's prime meridian at J2000. This yields: 





= R3{^^)Ri{-v)R3{-^rj)R3{-e)\ y \ , (12) 

inertial \ ^ / fixed 

where ipr^ is the planetary longitude of Mercury's vernal point at the observation 
epoch, i?3(— 0) the rotation matrix with respect to the z-axis, and Q the planetary 
rotation angle defined in eq. ljlOp . This yields the expression of a point on Mercury's 
surface in what is called here the Hermean Mean Ecliptic of J2000 frame which de- 
pends on the planetary rotation (rate) and obliquity, which are our main parameters 
of interest. The Hermean Mean Ecliptic frame (HME) is the frame in which we ex- 
press our observable and perform all our calculations. It is the hermean analog to the 
Earth's Mean Ecliptic of J2000 frame. 

In our simulation images are identified through their centre point expressed in the 
Hermean Mean Ecliptic frame. This does not mean that it is assumed that each centre 
point corresponds to a landmark. Instead this means that a landmark present in the 
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image can be precisely positioned with respect to the image centre. This choice has 
no influence on the results of the rotation experiment analysis, as in reality images 
will be taken with existing knowledge of available targets. Not every surface point 
can be used as a landmark. The image correlation process that will be used to match 
images of a specific surface region, needs features that are sufficiently well resolved, like 
steep crater walls or albedo spots with very clearly defined borders. Additionally the 
solar incidence angle on the surface plays an important part in guaranteeing sufficient 
contrast and limited shadowing in order to successfully match two images. The reasons 
behind this will be discussed in section I2.4.4i The images to be paired in this way 
do not need to be centred on the same surface spot, but only to have a sufficient 
amount of observed surface in common, which in turn has to contain enough adequate 
features for the image correlation process. As no current images of Mercury's surface 
match the high resolution expected from the MPO High Resolution Imaging Channel 
(HRIC), we assume that adequate target regions are distributed evenly about the 
planetary surface. The important quantities for our simulation are the timing of 
image acquisition and the volume of image data used for the parameter estimation, as 
we will show below, and not the precise location on the planet of the simulated targets 
used. 

The above procedure reduces the process of generating image pairs to the determ- 
ination of position differences quite simple: 

Axt, t, = Rs i^v) Rl i-V) R3 i-^r,) [R3 (-e(tl)) Xt, - R3 (-e(t2)) Xt,] . (13) 

where the indices ti and t2 correspond to the epochs of the two images to be compared. 
This quantity represents the observable in our simulations. 

In order to realistically simulate the restitution of the rotation parameters from 
these simulated observables we apply observational errors on each individual measure- 
ment of the position of a surface point in inertial space. 

2.4 Error Models 

The following error sources are considered: 

• satellite attitude error (composed of the star-tracker pointing error and the lack 
of knowledge of the orientation of the camera with respect to the spacecraft bus) , 

• satellite position error, 

• image timing error, 

• image correlation error. 

The first three errors refer to the inaccuracy in the determination of the spacecraft's 
absolute state and orientation, the fourth is due to the image correlation process, 
mainly because of the difference in illumination, altitude and orientation between two 
photographs to be correlated. 

2.4.1 Attitude error 

The attitude error represents the misalignment of the camera with respect to the 
expected nadir pointing position due to an error in the attitude determination or 
control. It is composed of a star tracker component and a camera component. The 
error in the knowledge of the pointing of the star tracker with respect to a celestial 
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reference frame is considered to be better than 2.5 arcsec ( Milani et all 2001 ). The 



misahgnment of the camera with respect to the spacec raft optical bench (i.e. due to 



thermal stresses ) is considered to be below 2.5 arcsec (jMilani et al.l . l200ll ). 

Those errors create a shift of the surface image effectively captured by the camera 
with respect to the centre of the nominal nadir image. We model the error through 
a randomly generated shift vector that gives the direction and the modulus of the 
misalignment of the pointing direction with respect to the expected/nominal nadir 
position. For the purposes of our simulations we use a uniform random distribution 
C/[0,27r] for the direction of the attitude error, and a normal random distribution 
A^[0, cj^] to model its amplitude. The attitude error amplitude cjx has been set to a 
value of 5 arcsec in this study, corresponding to a worst case scenario in which the 
camera pointing error is at its maximum. This creates a cone of uncertainty centred on 
the nominal nadir (ground track) surface in which we randomly position the observed 
image centre. To quantify the impact of different camera pointing error levels on the 
libration and obliquity estimation process we perform simulations with errors ranging 
from 2.5 to 5 arcseconds. The results of this comparison are presented in section H31 

2.4.2 Position error 

The satellite position error is due to a shift of the satellite centre of mass with respect to 
the expected position and results in a surface image shift that contributes to the error 
on the estimation of the planetary orientation and rotation rate. As for the attitude 
error, we express it through a random error vector with a Gaussian distribution around 
the actual position of the spacecraft centre of mass. The variance of this error can be 
determined from the precision of the satellite position estimation process, which varies 
according to the availability of spacecraft tracking. The spacecraft position will be 
determined through the MORE radio-science experiment. We consider two different 
error magnitudes depending on the tracking possibilities. The first case corresponds 
to a very small spacecraft position error of the order of 20cm (rms) and is applied 
whenever the BepiColombo MPO spacecraft is in view of and tracked by a DSN station 
on Earth. The second case corresponds to an error larger by 2 orders of magnitude 
(~ 10m) an d represents a wors t case scenario w henever the spacecraft is not tracked 
from Earth (jMilani et al.l . I2OOII : lless et al.l . l2009l ^. 



We can conveniently divide the error in terms of along-track, cross-track and ra- 
dial components. The radial component of the spacecraft position will be a lot better 
constrained than the components inside the plane perpendicular to it (by about one 
order of magnitude typically for other space missions). Additionally a small error on 
the spacecraft altitude will be of limited impact to our purposes. It will not signific- 
antly alter the surface area observed by the camera, resulting in a very small pixel size 
variation only. A radial error of 10 meters made from an altitude of around 1500km 
results in a MPO HRIC pixel size variation of ~ 0.1mm and is thus negligible for our 
purpose^. By contrast an along-track or cross-track error will shift the camera field 
of view by the same amount on the planetary surface independently of the spacecraft 
altitude, for example a 10m error in the along-track/cross-track plane will result in 
the same shift on the surface position of the target observed. We can further sim- 
plify our reasoning because the MPO orbit is quasi polar and the orbital period of 
the MPO is much smaller than Mercury's rotation, meaning that the spacecraft track 



^The altitude of 1500km coincides with the nominal apoherm altitude of the BepiColombo MPO 
of 1508km. (see table 

^The BepiColombo HRIC has a physical sensors with 2048 x 2048 pixels, and a total field of view 
(FOV) of 1.47°. 



10 



is almost parallel to meridians. The along-track error corresponds almost exclusively 
to a shift of the image in latitude and the cross-track error to a longitudinal shift. 
The two error components are expected to be distributed normally, to be of the same 
order of magnitude and be independent of each other. The Earth-based tracking of 
the BepiColombo mission is expected to be performed by only one tracking station. 
We have considered the MPO visibility from the ESA Cebreros station. 

Therefore we apply a normally distributed longitudinal and latitudinal shift to 
each image centre with a magnitude of either 20cm or 10m depending on the visibility 
conditions of the spacecraft from the ESA Cebreros tracking station. In order to assess 
the impact of a lack of Earth-based tracking at the times where the surface is observed, 
a worst case with a constant spacecraft positioning error of 10m has been considered, 
and is discussed in sec l4.5[ 

2.4.3 Time tagging error 

Each image is tagged with its capture epoch. The time tagging error indicates a 
discrete error due to the wrong time clock reference on board the spacecraft. It is 
possible that the on-board clock rate changes with time by an unknown amount in 
the time between two images to be compared. This error step is random and with its 
value is associated a spacecraft position error vector, defined as the product between 
the velocity vector of the spacecraft per the step time. This means that the shifting 
is directed along track. 

As the ellipticity of the MPO orbit is low we can approximate the MPO periherm 
velocity by its mean velocity: 



where a is the semimajor axis of the MPO. Using GM^ = 22032.09 km^s ^ this yields 
a mean velocity of 2.54km/s for the nominal semimajor axis value of a = 3393.7km. 
Over 1 millisecond this corresponds to a displacement of 2.54 meters which can be 
almost fully attributed to the along track component. Therefore all errors can be 
neglected if the on-board clock can be synchronised to better than 1/10 of a millisecond 
(corresponding to below 25cm). This level of clock precision is easily achieved, and 
therefore this error source has been ignored in our simulations. 

2.4.4 Image correlation error 

The image correlation error is related to the process of correlation needed to match 
images. This error depends on the difference between the respective image illumination 
and resolution (altitude), as well as on the presence of suitable surface features in the 
observed region. As said previously, no high resolution photographic data presently 
exists of Mercury's surface. For the purpose of this simulation we therefore postulate 
that suitable surface features are evenly distributed on the planetary surface, and 
sufficiently abundant for the experiment to be possible. 

An important part of Mercury's surface is covered with albedo spots, small features 
with a very different albedo than their surroundings. As the albedo contrast tends 
to disappear for a very high solar incidence angle (Sun close to the zenith), the loss 
of albedo features and overall contrast render the image matching techniques all but 
useless. The same can be said of the usual topographic image matching. When 
correlating images of the planetary surface of Mercury we have two main characteristics 




(14) 
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of the surface to turn to: sharp topographic features (hke craters and faults), and 
albedo spots. The rims of craters or potential large faults are very useful to precisely 
align overlapping images with each other, given that they are finely enough resolved. 
The main limiting factor for this is the crater morphology itself, over which we do not 
have any control, and the crater illumination. If the Sun angle is too grazing (low 
incidence angle) or too zenithal (high incidence angle) we will either have very long 
shadows or completely absent shadows. Without shadow it is impossible to correctly 
po sition the features on the s urface. 

Jorda and Thomas! (|200d l simulated the pattern matching of albedo features and 



craters for a Mercury orbiter and concluded that sub-pixel accuracy can be achieved, 
if images obtained at phase angles < 5° and > 55° are excluded (for albedo spots) 
and if the difference between the phase angles of the images compared/matched are 
smaller than 35° (for craters). In our work we use their surface illumination criteria 
to select favourable targets and associated images, while guaranteeing pixel-size image 
matching errors at most. Therefore in this study we consider a standard 1-pixel error 
for each image pairing, scaled to the pixel size of the lowest resolution image. Based 
on the pixel size of the HRIC camera this represents an error of 2.5 as, or 5 to 18 
meters on the surface for MPO altitudes between 400 and 1500km. 

All these illumination constraints are to our knowledge not fixed precisely in a 
definitive way , but are to be understood as worst cases. A more detailed study by 
Vacanti et al.l (j2005l ) specifically tailored to the BepiColombo mis sion using lunar im- 
ages a s well as synthetic images reached similar conclusions as Jorda and Thorna^ 

A recent experimen t using image correlation to measure the libration of Epimeth- 
eus bv lTiscareno et al. I (i009) achieves sub-pixel accuracy around 0.5 pixels. No limit- 
ations on the surface illuminations are made, as it is the movement of the shape of the 
body that is tracked and not individual details at its surface. If such a sub-pixel level 
can be reached for the BepiColmombo rotation experiment, the level of uncertainty on 
the libration and obliquity estimation could be further reduced. The image correlation 
error has the same effect as adding an additional camera misalignment. Therefore a 
detailed study of the effect of different camera pointing errors (performed in sec. (|4.5|) ) 
can be used to quantify the effect of a reduced image correlation error. 

It should be noted that the worst case image correlation error level of 1 pixel is as 
large as the pointing errors best case. 



Illumination constraints 

The following nominal illumination constraints are used in this study: 

• Minimum solar elevation angle on the observed surface region: 35° , 

• Maximum solar elevation angle on the observed surface region: 85° , 

• Maximum difference in illumination angle between two images, to be able to 
compare them successfully: 35° . 

All images used have to satisfy these conditions in order to be paired and used in the 
libration and obliquity estimation process. 



2.5 Required precision level 

What level of error on the obliquity and libration amplitudes is acceptable if we are 
to reach the BepiColombo rotation experiment mission objectives? The targeted 
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precision for the and ratios stated in the BepiColombo mission goals (see 
Milani et all bonih is A(%^) < 10% and A(^) < 1%. By using eq.© and eq.© 
defined earlier we can not only relate the errors on the planetary parameters C20, 
C22 , T\ and 7 to the resulting estimation of and ^ but also express upper preci- 
sion limits on the estimation of the rotation and orientation parameters posed by the 
mission requirements. Transforming those equations to isolate r] and 7 we get: 
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Using these relations we derive expressions for the upper error limits on both obliquity 
and annual libration amplitude that guarantee a successful experiment. We consider 
Mercury's orbital eccentricity e and Mercury's orbital plane inclination with respect 
to the Laplace plane / as error-free known quantities. Since the errors on the gravita- 
tional spherical harmonic coefficient s C2n and C22 are expected to be determined with 



relative accuracy better than 10"^ (jless et al.l . bnOfll 'l. the dominating factor limiting 



the precise determination of both quantities is the knowledge on the obliquity and 
libration amplitude value. With respect to the errors on the obliquity and libration 
amplitudes, the latter can be expressed in terms of errors on C /ME? and Cm/C as: 



A7 




(17) 



(18) 



Rivoldini etall JioO^) have calculated interior structure models of Mercury for 
several plausible chemical compositions of the core and the mantle. For all models 
they have computed the associated C /MB? and Cm/C values and thereby constrained 
their plausible value to the following intervals: 



C C 
0.325 < ^ < 0.37 , 0.3 < < 0.6 . 

~ MR^ ~ ~ C ~ 



(19) 



These values togeth er combined with the current best estimates for 7 and r] by 
Margot et aD (|2007l ) {rj = 2.11 ± 0.1 arcmin,7 = 35.8 ± 2 arcsec), yield the 
following worst case estimate: 



At? < l.Oas , 

A7<3.1as. (20) 

Increasing the precision requirement on ^ to 5% instead of 10% increases the expec- 
ted precision level on the libration amplitude to 7 < 1.4 arcsec (keeping the maximum 
error constraint on the obliquity identical). 
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3 Surface observability 



In this section we describe the influence of mission design constraints on the surface 
observabihty (useful surface for the rotation experiment) expressed in terms of surface 
coverage. Which parts of the surface can be observed at a given mean anomaly of 
Mercury depends on the orbital characteristics of the MPO. We here use the orbital 
elements defined by the European Space Agency as the BepiColombo baseline scen- 
ario (see table [5]). This orbit with a polar configuration coupled to the low altitude 
represents a favourable case for the photographic coverage of the planetary surface 
and so for the rotation experiment. 



Parameter 


Symbol 


Value 


semimajor axis 


a 


3393.7 km 


periherm altitude 


a(l — e) — 
a(l + e)-R^ 


400 km 


apoherm altitude 


1508 km 


eccentricity 


e 


0.163244 


inclination 


i 


90.0° 


argument of periherm 


UJ 


16° 


longitude of ascending node 




67.7° 



Table 2: Initial BepiColombo MPO nominal orbi tal elements in the M ercury equatorial 
system on 21st July 2020, 0:00 UT as defined bv lOarcia et~ZI (|2007l '). R^ = 2439.7km 



This mission scenario proposes a mission launch date around July 21st 2014 res- 
ulting in an orbit insertion around Mercury on May 21st 2020 and reaching the final 
MPO orbit described above on July 21st 2020. Science operations are expected to 
begin after a one month commissioning phase around Aug 21st 2020. In this study we 
start the simulation of the orbit on the 21st of July but start the surface observations 
only on the 21st of August, to account for the effect of the small orbital evolution over 
the month of the commissioning phase. 

Previous studies have shown the importance of the spatial orientation of the MPO 
orbital plane for the libration and obliquity determination. This is due to the 3:2 
spin-orbit resonance of Mercury and the fact that the MPO polar orbit's plane will 
not precess. An orbital plane with /3 = 0°, where /3 represents the angle between 
the Sun- Mercury per i apsis line and the orbital plane at Me rcury perihelion, has been 
shown by Jehn et al. (2004) and Sanchez Ortiz et al. ( 20061 ) to be the most favourable 
configuration for the obliquity and libration determination. Our simulations confirm 
this and therefore we only consider /3 = 0° in this study. 

Due to the proximity of Mercury to the Sun the BepiColombo MPO will be sub- 
jected to intense solar radiation. This poses a certain number of thermal issues, 
one of which in particular heavily impacts the spacecraft operations (Benkhoff, pers. 
comm.). During certain periods of Mercury's orbit around the Sun the solar radi- 
ation on the solar panels will cause them to overheat. At the basis of this lies the 
maximum allowable temperature of the solar panels, which cannot be crossed under 
penalty of permanent solar panel degradation. To limit the temperature of the panels 
their inclination with respect to the incoming solar radiation has to be modulated as 
a function of the Mercury true anomaly. Ironically it is when Mercury and the space- 
craft will be closest to the Sun, when the most energy is available for potential harvest, 
that the situation will be critical. The enhanced solar radiation available at perihelion 
could be insufficient to offset the loss of power due to the larger solar panel inclination 
than at further distances of Mercury to the Sun. The power available to the MPO 
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will then be set equal to the bare minimum needed to keep the MPO alive (thermal 
management, attitude control, minimal communications), and science operations will 
not be possible. We have implemented these science blackout periods for the MPO as 
additional constraints in our simulations. 

The primary science blackout windows considered are situated between 20° to 50° 
and 310° to 340° of Mercury true anomaly. These windows are symmetrically disposed 
around the Mercury perihelion, due to the fact that the nominal MPO orbit has its 
orbital plane aligned with the planet-Sun direction at Mercury perihelion (/3 = 0°). 
We study the presence of this exclusion zone as a worst case scenario, so as to predict 
its impact on the quality of the rotation and orientation estimation (librations & 
obliquity) in the frame of the BepiColombo rotation experiment. 

In the current ESA outline for the MPO science operations it is not planned to 
make use of the SIMBIO-SYS HRIC in the first 180 mission days. The HRIC is crucial 
to the rotation experiment, as it is the only on-board camera with a pixel size smaller 
than the libration and obliquity signals we strive to estimate. The HRIC camera 
consists of a 2048x2048 pixels detector with a field of view (FOV) of 1.47°, implying a 
pixel size of 5 meters at an altitude of 400 km and 19 meters from 1500 km. In order 
to evaluate the impact of a 180 day delay in HRIC operations on the feasibility of the 
Rotation Experiment we will also study short mission durations of only 180 days (the 
nominal mission being 360 days) . This will be discussed in section |H 



3.1 Nominal orbit 

In order to be able to determine libration and obliquity from camera images, land- 
marks of the planetary surface need to be flown over at least twice under favourable 
observation conditions. Therefore we determine the consequences of the various ob- 
servation constraints imposed by the rotation experiment on the maximum number of 
image-pairs obtainable over the whole planetary surface after the nominal mission dur- 
ation of 360 days and for the nominal orbit. Figure (jlap shows how many image-pairs 
can be generated after 360 mission days (nominal) for the nominal ESA orbit under 
the previously deflned nominal illumination constraints for each image-sized surface 
region on Mercury's surface. 

Observing at day-side is not a sufficient condition for having useful images. Because 
of the restrictions on the solar incidence (see section I2.4.4p observations at latitudes 
North and South higher than 55° where the Sun never rises above 35° in the sky are 
excluded for the determination of Mercury's rotation. We note that due to Mercury's 
very small obliquity the solar illumination angles are almost symmetrical with respect 
to the equator. The lack of image pairs at high latitudes is clearly visible in fig. p^ . 
Small patches of surface where the Sun is close to a zenithal position on overfiight 
will also be excluded due to the lack of contrast arising from this situation. This 
corresponds to the illumination condition imposing that the Sun be lower than 85 deg 
in the local sky. The longitudinal distribution of the number of image pairs on the 
planetary surface as seen in figure (jlap can be partly explained by Mercury's 3:2 spin- 
orbit resonance and the choice of a (quasi) polar MPO orbit. Taking into account 
that the plane of the polar orbit of BepiColombo's MPO will hardly precess, there 
are 12 possible overfiights of a given meridian in the nominal mission duration of one 
terrestrial year (~ 2 Mercury days). This corresponds to one target overfiight every 
half Mercury rotation or 29 days. At the times corresponding to these intersections, we 
need to check the illumination conditions. Useful images can only be obtained when 
the illumination conditions on the surface are fa vourable. Since Mercury's obliquity 
is close to zero (2.11 arcmin. iMargot et al. (|2003)) the illumination conditions can be 
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considered as controlled by the difference between the planetary rotation angle and the 
true anomaly. There will be (on average over a range of longitudes) 6 opportunities 
when the surface will be on the day-side (the 6 others happening at night). Because 
of the 3:2 spin-orbit resonance, the orbital configuration of Mercury repeats exactly 
every Mercury day (176 days). Therefore, for each observed target on the surface 
there are two groups of (on average) 3 images, which are identically illuminated and 
observed at the same Mercury true/mean anomaly. Since the libration in longitude 
is a function of the planetary mean anomaly, correlating images taken exactly one 
Mercury day apart yields no information on the annual libration amplitude 7. 

The image pair maps in fig. (|lap do not contain image pairs made at an identical 
planetary anomaly, nor images with illumination conditions differing by more than 
35deg of solar incidence, hence most of the correctly illuminated planetary surface is 
observed less than 6 times yielding between 1 and 5 image-pairs. Certain longitudes 
will be oversampled with respect to others because mission durations are not integer 
multiples of one Mercury day. 

Although the above reasoning is made for the case of a Keplerian orbit without any 
secular or periodic variations the image repeatability is changed only slightly by in- 
serting the perturbative effects of the low-degree gravity field harmonics. Nonetheless, 
as we have mentioned earlier a Keplerian orbit is not a good approximation since the 
perturbing effects of the low-order gravity field harmonics will be quite large relative 
to the signal of the rotation experiment parameters. Therefore, only on average after 
176 days the MPO will fly over the same surface regions and the whole set of possible 
observations will take 176 days or longer to fill (in steps of half Mercury rotation 
periods, corresponding to ~ 29 days, because the orbit plane does not precess and the 
planetary rotation unfolds under the satellite). 

In contrast to previous studies we do not consider a camera offset capability in 
order to in crease the number of observ ations of a given target / surface region. Sim- 
ulations by Sanchez Ortiz et al. ( 20061 ) have shown that a small camera offset capab- 
ility yields more precise libration and obliquity estimates for an identical number of 
observations, because more favourable target observations and otherwise impossible 
image-pairs can be obtained. At the present time the only camera offset possible for 
the BepiColombo MPO is through off-nadir pointing manoeuvres. Their amplitude is 
limited as the MPO's payload risks overheating if not closely aligned with the space- 
craft nadir. It should be noted that the addition of such a capability adds additional 
uncertainty on the spacecraft pointing direction and therefore the advantage of more 
image-pairs is not entirely expected to benefit the libration and obliquity estimation, 
unless it is completely impossible to obtain satisfying image pairs without an off-nadir 
pointing . As there are currently no high resolution maps of Mercury's surface avail- 
able this cannot be excluded, but seems highly unlikely. The results we present here 
represent a worst case scenario in which no off-nadir pointing is performed. 

The fact that the image-pair opportunities exhibit a latitudinal symmetry (North- 
ern and Southern hemisphere a both equally covered) (see fig fTa|) is due to the fact that 
the initial orbital configuration is optimised to assure equal coverage of the Northern 
and Southern hemispheres over the nominal mission lifetime taking into account the 
secular drift of the periherm. In the nominal ESA (ESAS) case the periherm will 
advance by about 33° per year, and starting at — 16°South will reach 16°North after 
one year. By advance of the periherm we mean the fact that the periherm occurs 
earlier in the orbit path than it would have been expected. 
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(a) Distribution of surface regions which are observable at least twice and from which at 
least one image comparison can be constructed. The colour scale represents the number of 
such image comparisons which are possible. The observations have been filtered according 
to an illumination constraint excluding all low solar elevations (below 35°), very high solar 
elevations (above 85°) and imposing that inside an image-pair the illumination angle 
difference cannot be greater than 35°. 
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(b) In addition to the same illumination constraints as fig. paf) . the observations have 
been filtered according to a science exclusion zone between 20° to 50° and 3f0° to 340° 
of Mercury true anomaly. 

Figure 1 

3.2 True Anomaly limitations 

Because of reasons of thermal management of the solar panels no science operations 
can be performed in two zones in true anomaly close to Mercury's perihelion (see 
above). Fig. ljlbp shows the effect of the addition of science exclusion zones between 
20° to 50° and 310° to 340° of Mercury true anomaly, on the repeated surface observ- 
ability. Comparing fig. fjlbp with fig. ljlap we notice two additional empty longitudinal 
bands situated between 20° and 50° of East and West longitude. The excluded true 
anomalies are directly related to the planetary longitudes affected due to the 3:2 spin- 
orbit resonance. Therefore we can say that if an exclusion zone contains the Mercury 
perihelion, the observations of the surface region centred on the prime meridian will 
be reduced or even completely vanish. The main effect of such a science blackout 
zone is to further reduce the available surface in which to look for suitable image-pair 
candidates. In the following section we compare the effects of the presence of such a 
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limitations on the libration and obliquity estimation process. 



4 Results for librations and obliquity 

In this section we investigate the quality of the estimation of the studied parameters 
(libration amplitudes and obliquity) for the nominal BepiColombo MPO orbit in order 
to establish the feasibility of the rotation experiment at the desired precision level. 
The following parameters are considered as limiting factors on the surface observab- 
ility: the solar incidence on the surface seen by the camera and the solar incidence 
variation between two measurements of the same target, the mission duration, the 
number of observed targets and the presence and size of science blackout zones. We 
studied the effect of various combinations of these parameters on the estimation qual- 
ity of the obliquity and libration amplitude. We also studied the effect of limiting the 
observations made by the MPO above 1000 km because of the trade-off bet ween the 
resolution and the amount of pictures ( Milani et al. . 2001 : Jehn et al. . 20041 ) . 



The results discussed in this section are obtained using simulated targets located 
inside the observable zones defined in the previous section. Using these filtered obser- 
vations we determine the value of the observable (the position change of a given image 
pair), add noise to the observable according to the error models and calculate the 
partial derivatives of the observable defined in eq. (jl3p with respect to the four para- 
meters studied here. This is followed by the estimation of the parameters through a 
least squares approach, solving the following system: 



p = A^C,-i6(A^C,-iA)-i , (21) 
-^u = -Q^, for i = l.A and / = ??, 7, Tif, 79 , (22) 

where A is a (m x n) matrix of the partial derivatives of our observable with respect to 
the different parameters (with n the number of parameters (n=4) and m the number of 
observations), A"^ is the transpose of A, 6 is a {m x 1) vector containing the residuals 
of the observations, and p {n x 1) is the common value taken by the parameters 
throughout all observations which is to be determined by our least squares approach. 
Ch is a (m x m) weight matrix that can be applied to the observations. Generally it is 
a diagonal matrix composed of the inverse a priori uncertainties on the observations 
(6) is used. In our case C;, is equal to the identity matrix as we do not give more 
weight to certain observations relative to others. A"^C^^ A is a square (n x n) matrix, 
and is called the normal matrix. 

We present here the formal errors reached on each parameter for different sets 
of constraints as a function of either the number of surveyed targets or the mission 
duration. 



4.1 Number of surface targets surveyed 

Figure ([5]) presents the formal error on each of the four parameters as a function of 
the mission duration for different numbers of surveyed targets (10, 25, 50, 100). The 
presence of a slow and steady decrease in formal error with mission duration for any 
constant number of surveyed surface targets is a strong indication that the data volume 
is not the bottleneck of the experiment. Indeed throughout the mission the surveyed 
targets will continue to yield additional image-pairs. The small decrease in formal 
errors when the target numbers increase is mainly due to the increase in data volume 
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and the fact that the formal error is proportional to where is the number of 
image-pairs. With a smaller number of targets the error on the rotation parameters 
increases but a sufficient precision could be obtained for a longer observation time. 
Nevertheless a large number of targets may be needed to reduce the error levels below 
the mission requirements. 

Results with the required precision (see eql20p can be reached after 180 mission 
days already for the obliquity and annual libration amplitude, for 50 surveyed targets, 
as all different observation windows will already have been explored after one Mercury 
day. For a small amplitude and long period signal such as both long period planetary 
librations studied here, the observation time series has to be long enough to cover a 
significant part of the cycle so that the cycle can be determined accurately. The main 
error sources of the experiment will have to be reduced in order to reach the required 
precision on the planetary libration amplitudes, as more measurements and a longer 
mission duration do not suffice (see figiH). A limited number of 25 targets suffices to 
reach an error level satisfying eq.(|20p on both the annual libration and obliquity in 
the nominal 360 day mission. 

All parameters except the planetary induced libration amplitudes can be solved 
for with a formal error smaller than the minimum mission requirement with 25 image- 
pairs after around 300 mission days as can be seen in fig.(j3|). The error on the Jupiter 
induced libration estimate is close to the mission requirements defined by eq. (pO]l after 
about 550 days. A doubling of the number of surveyed targets would enable to achieve 
a halving of the error value to barely reach the minimum requirement. 

4.2 Altitude limit 

If measurements are restricted to 1000 km in altitude (reducing the data volume) all 
the parameters are less well solved for, except for the libration amplitude (figl3|). The 
latter is best determined if observations are restricted to lower altitudes. This differ- 
ence is due to a combination of the following two facts: The annual libration signal is 
strongest in the equatorial regions as it is a longitudinal signal, and the periherm of the 
MPO orbit is located close to the equator creating the possibility of very strong-signal 
image pairs. If we allow images taken at higher altitudes we add image pairs of lower 
resolution spread over a wider latitude domain. Since we have not weighted the ob- 
servations as a function of the latitude or altitude or any other characteristic (because 
this would have skewed the results of the other three parameters while favouring this 
one), the addition of lower resolution image pairs degrades the total solution, diluting 
the signal coming from high signal to noise image-pairs taken at low altitudes. The 
solution therefore gets noisier. We don't see this for the other two libration amplitudes 
because of their very long period and therefore significantly smaller longitudinal dis- 
placement with respect to the annual libration movement. Their signal is weak enough 
to benefit from the presence of higher altitude observations, which enable a greater 
number of image pairs. The fact that the obliquity does not benefit from the low 
altitude only approach is explained in the following way. Due to Mercury's 3:2 spin 
resonance and the fact that the MPO orbital plane stays fixed in space and is initially 
fixed to (3 = Odeg, two of the three favourable dayside surface observations (see section 
[3]) happen symmetrically around the Mercury aphelion, and correspond to low alti- 
tude images taken from the MPO periherm orbit side. The third observation window 
is situated around Mercury perihelion and corresponds to high altitude observations 
from the MPO apoherm side. Because both low altitude observations happen almost 
symmetrically about perihelion their comparison will yield no information on Mer- 
cury's obliquity as the planet will be almost identically oriented in space with respect 
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to the MPO orbit. The third observation window provides a different viewing angle 
on the planet, allowing for image-pairs containing a strong obliquity signal. Therefore 
it is critical for the obliquity measurements not to limit the surface imaging to low 
spacecraft altitude. 

4.3 Science blackout zones 

In fig. ([3]) we compare the results obtained for a fixed number of 25 surveyed targets as 
a function of the mission duration, for different blackout zones and a case without any 
blackout zone. The zones considered correspond to the nominal case described earlier 
(see section [3|), one variant shifted towards Mercury perihelion, and one radical variant 
excluding all observations made between 300 deg and 60 deg of Mercury true anomaly. 
For this last zone we can see in fig.dH) some similarities with the altitude limitation as 
no observations can be made at and close to Mercury perihelion, such as no success- 
ful obliquity estimation. The results for all blackout zones allowing observations at 
perihelion are extremely similar and close to the result of unaffected case. This shows 
us that the size of these exclusion zone is not the dominant factor in the degradation 
of the parameter estimation process, but the presence of observations made close to 
perihelion is. Surface observations between —10° and -|-10° of Mercury true anomaly 
seem to be very important for the achievement of the mission requirements as the 10° 
to 40° blackout zone is barely able to achieve the necessary precision in the nominal 
mission duration of 360 days. 

4.4 Illumination constraints 

Figure [5] shows the formal error as a function of the mission duration for a fixed number 
of 25 observed targets for the nominal illumination constraints versus relaxed and 
tightened illumination constraints. The first goal of this figure is to show the absence of 
influence of the difference in illumination angles between images forming an image pair 
(separation limit) on the parameter estimation. The solid and dashed lines correspond 
to the nominal illumination conditions and the case of slightly relaxed conditions 
without the need for a maximum image illumination separation of 35 deg inside an 
image pair. Comparing them we cannot see significantly smaller errors for the dashed 
line. Comparing the nominal conditions to slightly stricter conditions represented 
by the dash-dotted line it can be seen that the stricter conditions enhance the errors 
significantly only for the libration amplitude. Finally the dotted line represents a quasi 
unlimited illumination scenario where only local solar elevations below 10° are being 
rejected. This allows for a hugely increased portion of the planetary surface, between 
the northern and southern 80° parallels and with an almost homogeneous longitudinal 
distribution, in which targets can be observed. Nevertheless for an identical number 
of surveyed targets, the nominal illumination conditions yield only about two to three 
times less favourable error bars on the parameter estimates. Again the annual libration 
amplitude is the most affected. We can clearly see that in all cases the formal errors 
decrease with mission duration. A satisfying precision level below 1 arcsecond for the 
obliquity and 3 arcsecond for the annual libration is reached after ~300 mission days 
for all sets of illumination constraints. From the perspective of the annual libration 
amplitude enough information is available to reach a precision smaller than defined by 
eq. ()20p after approximately 250 mission days, for any of the illumination conditions 
studied. 
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4.5 Influence of the camera pointing error level 

In order to gauge the impact of the camera pointing error level on the libration and 
obliquity estimation, we have simulated the experiment for error levels of 2.5, 3.8 and 
5.0 arcseconds. The comparison can be found in figure ^ and shows that a reduction 
of the pointing error level by a factor of two reduces the error level on the obliquity 
and libration amplitudes by slightly less than a factor of two. If the pointing error 
can be reduced to its minimum level of 2.5as the Jupiter induced libration could be 
determined within the nominal mission duration (360 days) with an error level below 
the requirements to pose an additional constraint on the planetary interior (< 2.5as). 
In general it can be said that a reduction of pointing error level from its worst case 
to its ideal case does not impact the order of magnitude with which the obliquity and 
libration amplitudes will be determined. 

The impact of the MPO position error on the rotation experiment has been studied. 
The nominal scenario in which the MPO is tracked by one Earth-based groundstation 
has been compared to a situation with a constant 10m position error. When tracked 
the MPO position error drops to 20cm. No impact on the final result from a constant 
10m error on the MPO position can be seen in fig.dU]). 

A position error of 10m on the MPO position translates in a maximum error of 
10m on the planetary surface. By contrast a pointing error of 2.5as from an altitude 
of 400km corresponds to a surface displacement of 4.8m. For S.Oas this corresponds to 
9.6m. From an altitude of 1500km this represents an error of 18m or 36m respectively. 

Two images composing an image pair are rarely both taken while the MPO can 
be tracked from Earth. Therefore the impact of a constant 10m tracking error on the 
final result is limited. This is especially true for scenarios with larger MPO pointing 
error levels. 



5 Conclusion &; Perspectives 

In the frame of this work, we have simulated the BepiColombo rotation experiment 
a nd its potential out come. The BepiColombo orbit characteristics proposed by ESA 



( (jGarcia et al.l . 120071 ) (see tabi2])) were used in order to show their impact on the 
quality of the parameter restitution using the libration experiment observables. By 
projecting the spacecraft orbit on the surface of the planet and considering surface 
illumination conditions required for the image matching process, the distribution of 
potential images on the surface has been obtained. The uncertainty on the satellite 
position and orientation (imposed by the precision of the star tracker and radio sci- 
ence observations respectively) is included by subjecting each potential image to a 
positioning error. In practise, the location and time of possible repeated photographic 
measurements of selected target positions on the surface of Mercury have been simu- 
lated and used as input for the determination of the rotation parameters. Additional 
errors arise when comparing two images of the same target region on the planetary sur- 
face. This process is limited by the photographic resolution of the least well-resolved 
image and depends on the solar illumination angle of the planetary surface region 
observed. 

The location and time of possible repeated photographic measurements of selected 
target positions on the surface of Mercury have been simulated and used as input for 
the determination of the rotation parameters. By means of a least squares process, 
three libration amplitudes and the obliquity are estimated. 

The accuracy of the reconstruction of the rotational motion of the planet has been 
studied as a function of the quantity of measurements made, the number of different 
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targets considered, the mission duration, the illumination conditions, the existence of 
observational altitude limitations and the size of potential spacecraft science blackout 
zones. Several parameters influence the quality of the libration and obliquity restitu- 
tion directly and indirectly through the limitation of the surface observability in time 
and space. Our results can be used to determine the sweet spot in observations and 
target numbers needed to achieve a good result while keeping the data volume to a 
minimum (because dedicated imaging opportunities will be limited), and can serve as 
a measure of the cost of additional observation in contrast to the gains obtained. 

As shown in section [3] the nominal image illumination constraints used in this 
study greatly reduce the surface regions available for target selection as well as the 
number of image pairs that can be made from these targets over a nominal mission 
duration of 360 days. This explains the relative insensitivity of the observable surface 
regions to changes in blackout-zone size, as their consequences manifest themselves 
mainly outside of the observable regions. In any case most image pairs are to be 
found at latitudes between 30° North and 30° South and in two windows centred 
around the planetary longitude 0° and 180°, due to the strict illumination conditions. 
Fortunately as the libration is a longitudinal motion, shifts between images are larger 
at lower latitudes. Therefore, images of near-equatorial regions will be of greater 
importance to the rotational experiment. 

The precision on the estimated libration amplitude also depends on the number of 
possible images of the same target used. More surveyed targets provide a larger data 
volume and more opportunities for high-signal image-pairs, but are limited by the data 
volume that will be allocatable to the rotation experiment from total BepiColombo 
mission data volume. 

The criteria for the best determination of the obliquity are different from those 
of the libration amplitudes because low altitude observations are near symmetrical 
with respect to Mercury's apohelion and therefore don't yield sufficient information 
on the obliquity in contrast with the annual libration. The incompatible optimal 
requirements for obliquity and libration set additional constraints on the strategy to 
obtain both quantities from this image based experiment, such as the requirement to 
use high altitude images taken on the MPO's apoherm side in order to be able to solve 
for the obliquity. 

We have shown that, in order to achieve the rotation experiment goals, if a science 
blackout zone has to be implemented it cannot encompass the planetary perihelion, 
under penalty of rendering the obliquity estimation impossible, and heavily degrad- 
ing all libration estimates. Surface observations between —10° and +10° of Mercury 
true anomaly are critical to correctly estimate Mercury's obliquity and in turn better 
constrain the annual and planetary librations. 

The two long period planetary induced librations can only be reliably solved for, 
if the annual libration amplitude and the obliquity are well determined. 

The quality of the spacecraft positioning is very important since it greatly in- 
fluences the ability to position the images taken by the spacecraft on the planetary 
surface and hence to estimate the planetary rotation. It is worth to remember that 
we have considered a worst case scenario for the magnitudes of each individual obser- 
vation error. This was done in order to guarantee not to be overly optimistic in terms 
of the quality of the estimation process because of the relatively basic error models 
used. 

For the nominal mission scenario we have shown that: 

1. The obliquity can be accurately solved for if surface observations around Mercury 
perihelion can be guaranteed. This is true even for a small number of 25 targets 
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surveyed. 



2. The annual libration amplitude is well solved for in almost all circumstances 
reaching error levels of around 1 arcsecond in true and formal error after less 
than 360 days and for data sets limited to 25 targets only. 

3. The main Jupiter and Venus induced libration amplitudes are best estimated for 
longer mission durations because of their long periods. 

4. The amplitude of the Jupiter induced libration is expected to be of the order of 
the annual libration amplitude, because the free libration period of Mercury is 
expected very close to the Jupiter induced libration period ( Dufev et al. . 20081 ). 



A good estimation of the Jupiter libration amplitude is possible under nominal 
conditions (as long as the obliquity is well resolved) and could deliver a constraint 
on the value of the free libration period. If the pointing error can be reduced 
to its minimum level of 2.5as the Jupiter amplitude could be determined within 
the nominal mission duration (360 days) to below < 2.5as. 

5. The Venus libration amplitude is ten times smaller than the annual libration 
amplitude and cannot (even for longer mission durations) be solved for accur- 
ately enough to deliver an independent constraint on the interior parameters of 
Mercury. The error levels never dip below 50% of the libration amplitude. 

6. All these results can be obtained when considering the worst case camera point- 
ing error level of 5.0as 

For all these reasons we can conclude that the nominal mission scenario as currently 
defined by ESA, taking into account the possible need for science blackout zones, 
is able to solve for the obliquity and annual libration amplitude and almost for the 
Jupiter induced libration amplitude, given that 25 targets are observed repeatedly 
over the nominal mission of 360 days, and that surface observations close to Mercury 
perihelion are possible. 
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Errors in function of mission duration, for different fixed target numbers 
Sigma on Obliquity (126.6 as) 
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Figure 2: Prom top to bottom the formal error (Icr) on the four parameters we 
solve as a function of the mission duration. Each plot contains four lines representing 
each a different number of targets used to obtain the solution. The thick black line 
represents the error level required for a successful experiment. The full, dashed, 
dash-dotted and dotted curves represent 10, 25, 50 and 100 targets used 
respectively. 
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Errors in function of number of targets surveyed, for different altitude limitations 
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Figure 3: From top to bottom the formal error (Icj) on the four parameters we 
solve for a 360 day mission for the nominal ESAS orbit. Each plot contains two curves 
representing a different upper altitude limit on the observations used in the results. 
The thick horizontal black line represents the error level required for a successful 
experiment. The black curve represents the altitude unlimAted scenario where images 
are not filtered based on the altitude from which they were taken. The dashed line 
corresponds to an altitude limitation of 1000km on the observations used. 
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Errors in function of mission duration for different science exclusion zones 
Sigma on Obliquity (126.6 as) 
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Figure 4: From top to bottom the formal error (Ifx) on the fom parameters we 
solve for a 360 day mission. Each plot contains four lines representing each a different 
Mercury true anomaly exclusion zone constraining the observations used. The thick 
black line represents the error level required for a successful experiment. The black 
curve corresponds to the absence of a blackout zone. The dashed, dashdotted and 
dotted curves represent the different exclusion windows studied. For the obliquity, the 
dotted line is always above 2 arcsec. 
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Errors in function of mission duration for different solar incidence limitations 
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Figure 5: From top to bottom the formal error (Ifj) on the four parameters we solve 
for a 360 day mission. Each plot contains four lines representing each different sur- 
face illumination conditions constraining the observations used. The thick black line 
represents the error level required for a successful experiment. The black curve cor- 
responds to the nominal illumination conditions. The dashed curve represents slighly 
relaxed conditions, without the illumination separation limitation for image pairs. The 
dashdotted curve corresponds to slightly more strict illumination conditions. Finally 
the dotted curve represents very lax illumination conditions without any separation 
limitation. 
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Figure 6: From top to bottom the formal error (lu) on the four parameters we 
solve for a 360 day mission. Each plot contains four lines representing each different 
camera pointing (in arcsec) and spacecraft positioning (in m) error levels. 

The black curve corresponds to the worst case pointing error scenario of Sarcsec. 
The grey curve represents an intermediate scenario, with 3. Sarcsec of pointing error. 
The dotted curve corresponds to the most optimistic pointing error level of 2. Sarcsec. 
The dashed curves stand for cases where the spacecraft positioning error is fixed to 
10m (black) or 20cm (grey). Finally the dashdottcd curve represents a case without 
pointing error, and with a constant 20cm position error. 
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